Local spin spirals in the Neel phase of L^-^Sr^CuC^ 
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Experimental observations of lightly doped I^-^Sr^CuCU, x < 0.02, revealed remarkable mag- 
netic properties such as the incommensurate noncollinear ordering (additional to the Neel ordering) 
and a tremendous doping dependence of the uniform longitudinal susceptibility. We show that the 
spiral solution of the t-t'-t" -J model obtained by taking into account the Coulomb trapping of holes 
by Sr ions describes these puzzling data perfectly well. Our solution firstly explains why the in- 
commensurate structure is directed along the orthorhombic &-axis, and secondly allows a numerical 
calculation of the positions and shapes of the incommensurate neutron scattering peaks. Thirdly, 
we calculate the doping dependence of the spin-wave gap, and lastly, we study the longitudinal 
magnetic susceptibility and show that its doping dependence is due to the noncollinearity of the 
spin spiral. 

PACS numbers: 74.72.Dn, 75.10.Jm, 75.30.Fv 75.50.Ee 



I. INTRODUCTION 

The phase diagram of La2- a; Sr a ;Cu04 (LSCO) shows 
that the magnetic state changes drastically with Sr dop- 
ing. The three-dimensional antiferromagnetic (AF) Neel 
order identifiedi below 325K in the parent compound 
La2Cu04 (LCO) disappears at doping x ~ 0.02 and gives 
way to the so-called spin-glass phase which extends up 
toift! 0.055. In both, the Neel and the spin-glass phase, 
the system essentially behaves like an Anderson insulator 
and exhibits only hopping conductivity. Superconductiv- 
ity then sets in for doping x > 0.055 (see Ref. 0). One 
of the most intriguing properties of LSCO is the static 
incommensurate magnetic ordering at low temperature 
observed in elastic neutron scattering experiments. Such 
ordering is a generic feature of LSCO, because it has 
been observed in the superconducting phase 3 -, in the spin- 
glass phas o 4 i 5 i 6 , and in the Neel phasai The incommen- 
surate ordering manifests itself as an additional elastic 
scattering peak shifted with respect to the AF position: 
Q = Qaf + 5Q, where Qaf = (tt, ±7t), setting the lat- 
tice spacing a equal to one. According to experiments in 
the underdoped region 0.055 < x < 0.12, the shift scales 
linearly with doping and is directed along the crystal axes 
of the square lattice, 5Q « 2x(±ir, 0) or SQ 2x(0, ±7t), 
see Ref. U In the spin glass phase, 0.02 < x < 0.055, the 
shift also scales linearly with doping, but it is directed 
along t he o rthorhombic b- axis, (5Q w ±V2x{ir,-ir), see 
Refs. I ll lsltil Finally, in the Neel phase for x < 0.02, the 
incommensurability is almost doping independent and di- 
rected along the orthorhombic 6-axis (analogous to the 
spin glass phase), <5Q w ±0.02v2(7r, — tt), see Ref.0 Ex- 
perimental data for the elastic neutron scattering prob- 
ability at x = 0.01 and x = 0.014 are shown in Fig. 3 of 
Ref. □ (and reproduced in Fig. [8] in this paper). The ex- 
tracted correlation lengths summarized in Tab. 1 of Ref. 7 
clearly indicate the presence of long-range incommensu- 
rate correlations. The correlation length is about 200 A 



along the orthorhombic ^-direction and more than 500 A 
along the a-direction. To resolve the small incommensu- 
rate peaks in the background of the huge commensurate 
peak, the authors of Ref. used the fact that the scat- 
tering amplitude for neutrons interacting with electron 
spins is of the form 



T q = M W • B q = 4tt/x 



(IV) 



q (Mq • q) 



(1) 



/Xq, inconsistent 



Here is the magnetic moment of the neutron, q the 
momentum transfer and B q the Fourier transformation 
of the magnetic field generated by the magnetization den- 
sity \x (r), which in momentum space becomes /x q . It is 
well known 2 ", that due to the Dzyaloshinski-Moriya (DM) 
interaction, the commensurate Neel magnetization is di- 
rected along the orthorhombic 6-axis, (/x q ) oc (1, —1). If 
one choses the momentum transfer q along this direc- 
tion, as has been done in the experiment**, Eq. Q shows 
that in this case, the commensurate magnetization does 
not contribute to the scattering, which allows one to ob- 
serve the small incommensurate peaks. This also clearly 
indicates that the incommensurate peaks are due to a 
noncollinear spin structure, T q oc /i^ 
with any collinear spin stripe picture. 

The doping dependence of the DM spin-wave gap in 
the Neel state has been measured quite recently^. The 
observed reduction of the gap is clearly due to the loss of 
Neel order, which is completely destroyed at 2% doping. 

Another remarkable feature of LSCO is the doping de- 
pendence of the uniform magnetic susceptibility in the 
Neel state at zero temperature. According to Ref. 0, the 
longitudinal susceptibility Xb changes tremendously al- 
ready at x = 0.01, while the transverse susceptibilities 
Xa and Xc remain practically unchanged^ compared to 
the undoped compound. 

In the present work, based on the spiral solution of the 
extended t-J mode l 11 ! 12 , we explain and calculate the 
magnetic properties described above in the lightly doped 
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Neel state and show why the incommensurate structure 
in the insulating state (in both the Neel and the spin-glass 
phase) is directed along the orthorhombic 6-direction. 



The idea of spin spiral formation in an antiferromag- 
net with itinerant particles goes back to Nagaokai^ who 
noted that for a sufficiently small superexchange a mobile 
particle turns the antiferromagnet into a ferromagnet. In 
relation to the t-J model, this idea was first formulated 
by Shraiman and SiggiaM, who pointed out that for an 
appreciable superexchange, it is energetically favorable to 
allow the collinear Neel state to relax and form a spiral, 
in which holes can hop more easily. For a long time, the 
issue of stability of the spiral state remained rather con- 
troversial, because of the negative compressibility of the 
uniformly doped statei£*i§iiLi2ii2i2ii. It has been demon- 
strated recentlj*ii that the next-nearest neighbor hopping 
matrix elements t' and t" are crucially important for the 
stability of the uniformly doped state and that the spi- 
ral is indeed stable for physical values of t and t'. In 
the uniformly doped state, the spiral is always directed 
along the crystal axes of the square lattice. The possi- 
bility of spiral ordering in the insulating spin-glass phase 
of LSCO has been pointed out in Ref. |2l|. In the insu- 
lating phase, x < 0.055, the compressibility issue is not 
important because of the trapping of holes by Sr ions. 
Since trapped holes induce a spiral directed along the 
diagonal of the square lattice^, the direction of the in- 
commensurate structure is rotated by 45° at the point of 
the percolation-like insulator-superconductor transition. 
In addition, the anisotropy of the dc- conductivity in the 
spin-glass phased has been explained and calculated in 
Ref. El 



Our paper is organized as follows. We first recall 
the main results of the extended t-J model concerning 
the single-hole dispersion and the Sr-hole bound state in 
Sec. Inland explain why the incommensurate spin struc- 
ture (spin spiral) is directed along the orthorhombic in- 
direction. In Sec. IIIII we introduce the effective low- 
energy Hamiltonian for doped antiferromagnets within 
the framework of the non- linear tr-model (NLSM) , which 
is a very convenient technical tool to deal with spin de- 
grees of freedom in the t-t'-t"-J model, especially when 
taking into account the magnetic anisotropics due to the 
DM and XY interactions. Numerical simulations of this 
model, presented in Sec. II VI allow us to study long-range 
correlations in the ground state at zero temperature. A 
careful comparison between our calculations of the in- 
commensurate neutron scattering peaks and recent ex- 
perimental results is contained in Sec. E] The evolution 
of the DM induced spin-wave gap upon doping is pre- 
sented in Sec. I VII Sec. lVIIl is then devoted to the doping- 
dependencc of the uniform magnetic susceptibilities and 
finally, we present our conclusion in Sec. IVIIII 




(a) (b) 

FIG. 1: (a) Schematic drawings of a GuOi plane in the 
tetragonal and orthorhombic phase, (b) Reciprocal lattice with 
tetragonal and orthorhombic (for simplicity shown without 
distortion) unit vectors. Due to the orthorhombic distortion, 
holes around St = (§, — § ) have lower energy and the sptrals 
induced by these holes are directed along the orthorhombic b- 
direction. This leads to the incommensurate magnetic peaks 
shown as small open circles. 



II. THE Sr-HOLE BOUND STATE 
("IMPURITY") AND PINNING OF THE SPIRAL 
DIRECTION TO THE ORTHORHOMBIC b-AXIS 

Over a decade ago, the 2D t-J model has been sug- 
gested to describe the essential low-energy physics of 
high-T c cupratesSiSSiSS. In its extended version, this 
model includes additional hopping matrix elements t' and 
t" to next-nearest neighbors. The Hamiltonian of the 
model is well known, see e.g. Ref. [llj, and we do not 
present it here. The numerical values of the parameters 
of the t-t'-t"-J model corresponding to LSCO follow from 
Raman spectroscopy^! and ab-initio calculations 28 . It is 
convenient to measure all energies in units of </, i.e., we 
set J = 125 meV — ► 1 and obtain t — 3.1, t! = —0.5 
and t" — 0.3. At zero doping (no holes), the extended 
t-J model is equivalent to the Heisenberg model and de- 
scribes the Mott insulator LCO. Removal of a single elec- 
tron from this Mott insulator, or in other words injection 
of a hole, allows the charge carrier to propagate. Single- 
hole properties of the t-J model are well understood 29 . 
The main features are a very flat dispersion along the 
edges of the magnetic Brillouin zone (MBZ) with four 
degenerate half-pockets centered at S — (±-|,±-|). The 
quasi-particle residue at the minimum of the dispersion 
is Z as 0.3. In the full-pocket description, where two 
half-pockets are shifted by the AF vector Qaf, the two 
minima are located at S a = (S, |) and Sb = (f , — f ), 
see Fig. db). The system is thus similar to a two- valley 
semiconductor . 

Let us now consider a single hole trapped by the 
Coulomb potential of the Sr ion and refer to this bound 
state as the "impurity". Since such an "impurity" is an 
intrinsic part of LSCO, the word could be misleading 
and we therefore use quotation marks to avoid confu- 
sion. An "impurity" has a hydrogen-like structure. In 
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coordinate space, one part of the hole wave-function cor- 
responds to the Is bound state and depends smoothly 
on r, while a second part, dictated by the Bloch the- 
orem, is rapidly varying with r and different for holes 
located near S a or Sj,. A pseudo-spin finally indicates 
if the hole sits on the |- or j-sublattice. The ground 
state can therefore be described by two quantum num- 
bers, the hole-pocket and the pseudo-spin associated with 
the sublattice. In this case, the bound state in the tetrag- 
onal phase is four- fold degenerate {ls a |i lsaj.; ls&f> Is&j,}- 
However, in the low-temperature orthorhombic phase of 
LSCO we are interested in, the orthorhombic 6-direction 
is slightly longer than the a-direction, as illustrated in 
Fig. Ilia), a = 5.349 A, and b = 5.430 A, see e.g., Ref.0 
To first order, this distortion only influences the diago- 
nal hopping matrix element t', whose contribution to the 
dispersion is equal to 

e q = Zt' a ( e <9-+*fl« + H.c.) + Zt' b ( e -^+ l «» + H.c.) . 

Using the estimations t' a = t' oc 1/a and t' b oc 1/b, it is 
easy to see that the degeneracy associated with the pock- 
ets is lifted and holes in the vicinity of Sb have slightly 
lower energy 

Ae = e s - e Sb = iZ(t' b - t' a ) ~ 4Zt'(a/b - 1) - 1 meV . 

(2) 

The bound state is thus only doubly degenerate 
{lsbi, 1 S 6J.}- The spiral formed by the trapped hole 
("impurity") is orthogonal to the corresponding face of 
the MBZ 12 . Hence this immediately explains, why in 
the orthorhombic phase spirals are directed along the 
&-direction. This consideration is equally applicable to 
both the Neel and the spin-glass phase. The spiral pin- 
ning energy is given by E p i n ~ £ 2 xAe, x being the hole 
concentration and £ the correlation length in units of lat- 
tice spacings. The correlation length is anisotropic and 
temperature dependent, and due to frustration it remains 
finite, even at zero temperature, see Ref. Il2l Using the 
experimental valued at x ~ 0.02 obtained at T = 1.6K, 
£ 2 ~ 138, one finds E pm ~ 3 meV ~ 35K. We would like 
to emphasize that Eq. J5J) only provides a crude estimate, 
an accurate LDA calculation of t' b and t' a is necessary to 
determine the precise value of Ae. 

III. EFFECTIVE LOW ENERGY 
HAMILTONIAN OF THE DOPED 
ANTIFERROMAGNET. THE NLSM APPROACH 

The NLSM is a very convenient tool to describe the 
low energy dynamics of the weakly doped t-J model, rep- 
resenting underdoped cuprates. In this framework, the 
staggered component of the copper spins located in a sin- 
gle layer of LSCO is represented by a continuous vector 
field n(r) of unit length n(v) 2 = 1. To avoid confu- 
sion, we denote vectors acting in the three-dimensional 
(3D) spin space by arrows and vectors acting in the 2D 
coordinate space by the usual bold font. Throughout 



this paper, we adopt the orthorhombic coordinate sys- 
tem shown in Fig. ^ with unit vectors e a , a = a,b,c. 
Coordinate and spin space are linked through the pin- 
ning of the commensurate Neel magnetization to the or- 
thorhombic 6-axis. It is therefore convenient to use the 
same coordinate system in both cases with real space unit 
vectors e Q = e a . 

An elegant way to incorporate the DM and XY 
anisotropies in the NSLM, crucially important for an ac- 
curate description of the Neel state of LSCO, has been 
proposed quite recently by Silva Neto et at. In the static 
limit, the energy of the spin system reads^i 

En = ^tJ (W(r))2 

+ (^-n(r)) 2 +r c n 2 (r)| . (3) 

Here xj_ ~ 0.065 is the magn etic susceptibility of the 
Heisenberg model, c = a/ p s /x± ~ 1-66 the spin- wave ve- 
locity and p s s» 0.18 is the spin stiffness. The anisotropies 
are due to the DM interaction, with a DM vector directed 
along the orthorhombic a-direction D m 0.02e a , and the 
XY-term, which leads to \fT~ c w 0.04. 

The trapped hole has a hydrogen-like ground state 
wave- function 

ijj(v) = *x(r) = *y^e- Kr , 

where VP is a two-component spinor (independent of r) 
describing the pseudo-spin. Note that according to the 
previous section the hole resides in the Sb hole-pocket, 
however, we do not write this index explicitly in the wave 
function. The value of k slightly depends on doping. It 
decreases with increasing doping, because of the Coulomb 
screening. According to data on hopping conductivity^ 
K ~ 0.3 — 0.4 at x — 0.002 and very recent preliminary 
data taken at x = 0.03 indicated that k ~ 0.2. The hole 
energy inside the "impurity" can be written as 

£20 = J d 2 r ^ (r) (-^ ~^Ai> to , (4) 

where (3 w 2 is the inverse mass of the trapped holsii, and 
e e is the effective dielectric constant. Strictly speaking, 
one has to write the Coulomb energy as e 2 / (e e V d 2 + r 2 ), 
where d ss 2 Ais the distance from the Cu02 plane to the 
Sr ion, but since r ^S> d one can safely neglect the d 2 
term. 

The interaction energy of a hole with the spin degrees 
of freedom readsi^ 

= V2g J d 2 r [n (r) x ^ (r) aip (r)] (e ■ V) n (r) . 

(5) 

Here e is a unit vector orthogonal to the face of the MBZ, 
e = e&, since the hole resides in the Sb hole-pocket and a 
are the Pauli matrices. The coupling constant g w Zt ss 1 
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has been calculated previousljiii*i£ within the t-t'-t"-J 
model. 

The restriction to the lightly doped Neel phase allows 
us to expand the n-field around the dominating commen- 
surate order directed along the 6-axis 



n (r) = e b + 7? (r) = e b + ^ n a (r) e c 



(6) 



where 7? is a small vector field orthogonal to the b- 
direction, i. e. 7? <C 1. Using this approximation, the 
total energy of a system with N "impurities" can be ex- 
pressed as 



E = E n + E^ + E,j, 
d 2 r 



(7) 



/ 



jf E {[V7r«(r)] 2 + M^[7r a (r)] 2 } 

^ ct—a.c 
N <| 

+v / 25E'°( r - r i)»"i (e- V)vr(r) ^+B , 

i=l J 

where m, = e), x f^ja^j^ are the directors of the 
"impurities", \rhi\ = 1, p(r) — x 2 (r), and Eq = 

N 2ne 2 ^j . Since Eg is independent of m 8 and 7?, it 

only leads to a shift in energy which can be safely omit- 
ted for the purpose of the present work. The "masses" 
follow from (J3J): 



M„ 



D/c = 0.012 and M c 



r c / c = 0.024 



Note that this definition of the "masses" differs from the 
standard one by a factor 1/c. 

A. Solution for an isolated "impurity" 

Let us for clarity reasons introduce a special notation 
for the single "impurity" problem and use the vector field 
9 instead of 7? in this case, i.e., analogous to the defini- 
tion JBJ, we use 



n (r) = e b + 9 (r) = e b + 



a— a,c 



and assume « 1. According to Ref. {H, the single 
"impurity" solution which minimizes the total energy Q 
satisfies the equation 

(-A + Ml) 9 a (r) - m a M(e • V)p(r) = . (8) 

Here m a is the a-component of the director m, and 
M. = \/2g/p s sa 8 is the effective dipole moment of the 
"impurity" . 

The solution of Eq. (JSJ in coordinate and momentum 
representation reads 

V (r) = m aJ —— {(1 + 2kt) e~ 2 ^ - M a rK x (M a r)} , 



2tt r 2 



(9) 



where K-y is a modified Bessel function of the second kind 
and we make use of the fact that M a <§; k. The Fourier 
transformation p q of p (r) is given by 



8k j 



(4K 2 + g 2 ) 3/2 



1 (q <C «) 



(10) 



The above solutions are valid in the case of an isolated 
copper-oxide plane. 

Influence of neighboring CVO2 planes 

Let us consider the modification of the solution J3J due 
to the interaction with other copper-oxide planes. The 
energy spectrum of spin waves in LCO reads^ 

w {q x ,q y ,q z ) = 2.32 J ^ (1 + a ± /2) 2 - (7,, + a ± j ± /2) 2 

7|| = 2 ' cos qx + cos OS'i ry± = cos 9z ' 

Here ~ 5 • 10~ 5 describes the superexchange between 
the planes. We disregard DM and XY anisotropies, be- 
cause they are already taken into account in the effective 
NLSM. Incorporating the q z dependence of this spectrum 
in Eq. @ in the limit where q = (q x , q y ) is small, we get 



a 



ql + q 2 y + 4a± sin^( gz /2) + M, 



72 



(11) 



In order to find an effective expression of (|ll|l in the 
single-layer approximation used to derive Eq. @, we in- 
tegrate ( flip over all momenta q z and obtain 



dq z 



* Jo g 2 + g 2 + 4a ± sin 2 (g z /2) + M2 



^(q 2 x +q 2 y + Ml) (4a ± + g2 + g2 + M 2) ' 

which can be approximated by the original expression 
with the substitution of an effective mass 



M a -> Mf = M a I 1 



4o_l\ 



1/4 



i.e. Mf = 0.015 and M c cff = 0.026. We use these values 
for further calculations. 



B. Solution for multiple "impurities" 

It is straightforward to generalize the single "impurity" 
solution to a system consisting of N "impurities" . 
Because of the assumption 7? 2 <C 1, the 7?-field is the 
linear superposition of N independent local spirals 



7? (r) 



JV 

cx.—a.c i—1 



(12) 
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where 9 a (r) is the single "impurity" solution given by 
Eq. ©. 

As already pointed out in R.ef. lT3. the single "impurity" 
solution is degenerate with respect to the orientation of 
the director m, since minimizing the total energy J7J only 
enforces the (9-field to be parallel to fa. This degener- 
acy is lifted by the interaction with other "impurities" . 
Having N holes trapped by N Sr ions, one can calcu- 
late the effective interaction energy between "impurities" 
due to perturbation of the NLSM in the limit where the 
distance r between two "impurities" is sufficiently large, 
r = |rj — Yj\ > 1/k ~ 3. Let us substitute l(T2|) in the 
total energy Q and split the result into an effective in- 
teraction energy U (containing contributions i ^ j) and 
a self-energy £ (arising from contributions i = j). Omit- 
ting terms independent of m, which just shift the energy, 
we find 

N N 

e = u+^= Y, u ij+ E ( 13 ) 

ot=a,c i^j a—a.c i 

where 

U % = E ^-M 2 a m?m«{(e l ■ e 3 )K (M a r) 
(ei ■ &j) 



2(e 4 • r)(ej • r) 



K 2 (M a r) } , (14) 



with the modified Bessel functions K n , and 



£? 



— M*(m?) In 



(15) 



For M a -> 0, Eq. O agrees with Refs. I2ll35ll3l We 
would like to stress that even though the interaction 1|14|) 
looks similar to the usual electrostatic dipole-dipole in- 
teraction, it is substantially different: Firstly, the dipole 
director rhi which determines the polarization of the spi- 
ral is decoupled from the vector ej which determines the 
direction of the spiral. Secondly, the sign of i fTljl is op- 
posite to what one naively expects from the analogy to 
electrostatics and thirdly, the interaction is of finite range 
due to nonzero masses. Although the self-energy (|15fl is 
proportional to and therefore small numerically, our 
results presented in Sec. show that its contribution 
is crucially important and leads to an alignment of the 
"impurity" directors rhi along the a-direction, because 
M a < M c . 



1. "Molecular impurities" 

The effective interaction l|14[) is valid for "impurities" 
separated by a distance r > 1/k. In this case, the or- 
thorhombic distortion favors holes to reside in the vicin- 
ity of Sb, and the vectors e; in l|14|) are all directed along 
the 6-axis, i.e., e$ = e^. However, as soon as two "im- 
purities" are sitting very close to each other, the holes 



can no longer be treated as independent particles but 
have to obey the Pauli exclusion principle, as explained 
in Ref. 0. In order to elucidate this mechanism in more 
detail, it is convenient to distinguish between "atomic" 
and "molecular impurities" : an "atomic impurity" is just 
a redefinition of the bound state formed by a Sr ion and 
the trapped hole (in analogy with the hydrogen atom), 
whereas a "molecule" describes two or more "atomic 
impurities" with noticeably overlapping hole wave func- 
tions. The formation of "molecules" is entirely due to the 
Pauli blocking: For two or more well separated "atomic 
impurities" , the Pauli exclusion principle does not apply 
and the orthorhombic distortion favors holes residing in 
the pocket centered at Sb- But as soon as "atomic impu- 
rity" wave-functions have a non-negligible overlap, the 
Pauli blocking sets in and two holes in the same hole- 
pocket must have opposite pseudo-spins, which prevents 
the formation of a local spiral and therefore does not 
lead to a gain in energy. In this situation, it will be en- 
ergetically favorable to place the holes in different pock- 
ets, which, due to the orthorhombic distortion costs an 
energy Ae « 1 meV J2J but allows the formation of a 
local spiral^. For parallel alignment of the "impurity" 
directors this gain in energy is of about 3 — 5 meV. A 
"molecule" can therefore be represented by the director 
m and the vector e = e a ± e;,, since it is a superposi- 
tion of two holes in pockets S a and Sb with two times 
the self-energy (|15|l of an "atomic impurity". The crit- 
ical distance between "atomic impurities" , below which 
"molecules" are formed, is not well defined. Clearly, it is 
smaller than the average distance between "impurities" 
at the insulator-superconductor transition (percolation) 
at doping x = 0.055, i.e. r c < \j\fx ~ 4.3 lattice spac- 
ings. Although the formation of "molecules" is crucially 
important to explain the jump of the incommensurability 
direction at the insulator-superconductor transition^, it 
is negligible in the lightly doped Neel phase considered 
in this work. For instance, for r c — 2 at doping x = 0.01 
only 6% of the "impurities" form "molecules" - a contri- 
bution which can be safely neglected. This conclusion is 
supported by extensive numerical simulations for r c < 3 
lattice spacings, which do not indicate any modifications 
due to "molecules" . 



IV. GROUND STATE AND CORRELATION 
FUNCTIONS 

Let us now investigate the properties of the zero- 
temperature ground state of the system. The basic as- 
sumption behind the derivation of the effective interac- 
tion (|13Jl is the existence of a dominating commensurate 
Neel order 10, which restricts our analysis to the low- 
doping situation x ^ 0.02. From a technical point of 
view, the Hamiltonian i|13|) represents a set of interact- 
ing dipoles m.j. In the static limit, where any kinetic 
energy terms are absent, we are dealing with a classical 
problem. This classical approximation is justified by the 



large value of the effective dipole moment M. w 8 and by 
the long-range character of the interaction. In order to 
describe realistic experimental situations, with samples 
of LSCO consisting of many copper-oxide layers, each of 
them containing a random arrangement of "impurities" 
and all of them contributing to the measurement, we have 
to average the quantities we calculate over many differ- 
ent realizations of random dipole positions. To find the 
ground state of the Hamiltonian (fH3|> . we perform clas- 
sical zero-temperature Monte-Carlo simulations with up 
to N = 200 randomly distributed dipoles on a square 
lattice. These dipoles are separated by an average dis- 
tance I = l/y/x lattice spacings. In accordance with the 
previous discussion of the spiral direction pinning, we 
set ej = e^. Our results clearly indicate that the self- 
energy term (|15f) leads to a pinning of the dipoles along 
the orthorhombic a-direction, TOj cx ±e a , because of the 
smaller mass. This effect allows us to divide the task of 
finding the ground state into two parts: in a first step, 
we only consider the Ising-like situation, where all dipoles 
are aligned along the a-axis. Starting from random initial 
conditions, we find the ground state of a given realization 
by exactly minimizing clusters of eight dipoles at a time. 
Our algorithm generates random walks through the sys- 
tem and forms the clusters at a given site according to the 
strength of the interaction l|14l) with neighboring dipoles. 
The total energy is then minimized with respect to the 
dipoles in this cluster and the algorithm proceeds to the 
next site of the walk. In a second step, we perturbate 
this collinear arrangement and allow the dipoles to relax 
and have a nonzero component along the c-direction. Wc 
find that for a given number of realizations, the percent- 
age of ground states with noncollinear dipole alignment 
decreases with increasing system sizes. However, this ob- 
servation could very well be influenced by the difficulty 
of finding the optimal noncollinear dipole arrangement in 
larger systems. 

A typical ground state configuration for a given real- 
ization of random dipole positions at x = 0.01 is shown 
in Fig. [3 One clearly identifies domains with paral- 
lel dipole alignments. The size of these domains along 
the a-direction is substantially larger than along the b- 
direction. This picture also reflects the ferromagnetic 
(resp. antiferromagnetic) character of the dipole-dipole 
interaction (|14(l in the a-direction (resp. o-direction). 
Having found the ground state of a given realization in 
terms of the "impurity" directors to, one can easily calcu- 
late the resulting n-field, using Eqs. ©, l|T2|) and As 
an example, Figs. Ola) and (b) show the n-field derived 
from the ground state dipole arrangement for dopings 
x = 0.005 and x — 0.014, respectively. For readability, 
we only show a small part of these systems, where "impu- 
rities" are represented by filled circles. One can see that 
the initially dominating orientation along the 6-direction 
becomes weaker with increasing doping. 

In order to characterize the ground state, we define a 
correlation function, closely related to the neutron scat- 
tering cross section, see Eq. (|17|) in the following Sec. 
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FIG. 2: Example of a ground state dipole arrangement for 
1% doping. The dipoles form clusters stretched along the or- 
thorhombic a-direction. 



D a = 
q 



iq-(r 



(16) 



Due to the very strong pinning of the dipoles along the a- 
direction, their c-components are negligible and the cor- 
responding correlations vanish, i.e., D c (q) = 0. The 
asymptotic behavior is thus given by D a (q) — > 1 for 
a — ► oo . The common feature of D", obtained for dif- 
ferent hole concentrations a; is a pronounced peak cen- 
tered on the orthorhombic 6-axis, as shown in Fig. Q] for 
x = 0.01. The statistical average of is obtained from 
200 realizations of random dipole distributions in sys- 
tems with L — 141 x 141 sites and N = 200 "impuri- 
ties" . Statistical fluctuations are negligible compared to 
these well pronounced incommensurate peaks, which are 
manifestations of the dipole domains shown in Fig. [21 A 
density plot of D a (q) (same data as Fig. shown in 
Fig. [3] reveals the anisotropy of the peaks. The width 
along the 6-direction is around twice as large as along 
the a-direction. Most likely, we slightly overestimate the 
width in the 6-direction and substantially overestimate 
the width in the a-direction due to finite-size effects. The 
substantial underestimation of the correlation length (in- 
versly proportional to the width of the peak) along the 
a-direction is clearly illustrated in Fig. The size of 
the domain along the a-direction is comparable to the 
size of the lattice. Unfortunately, the convergence of our 
Monte-Carlo procedure is getting very slow when simu- 
lating larger systems, so that we are limited to this size. 
In Fig. [21 we show slices through the maximum of the 
peaks parallel to the 6-axis. In this lightly doped regime, 
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FIG. 3: n-fields derived from two realizations of random impurity distributions at dopings x — 0.005 (a) andx — 0.014 (b). Both 
systems consist of L — 60 x 60 sites and contain N — 18 and N = 50 dipoles, respectively. For readability, only a characteristic 
part of these systems is shown, with circles denoting the positions of "impurities". Spin and coordinate space are linked by the 
pinning of the commensurate magnetization along the orthorhombic b-direction. Clearly, doped holes lead to a destruction of 
Neel order, experimentally observed at x = 0.02. 





FIG. 4: (Color online). Normalized correlation function 
obtained from ground states of model for systems with 
L = 141 x 141 sites and N = 200 (x = 0.01) "impurities". 
The statistical average is taken over 200 realizations. Due to 
the pinning along the a-direction, — 0. 



the correlation length £{, clearly decreases with increas- 
ing doping. Note that the correlations length extracted 
from neutron scattering experiments^ cannot be directly 
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FIG. 5: Density plot of shown in Fig. ^] The plot clearly 
reveals the asymmetry of the incommensurate peaks. The 
width along the b-direction is around twice as large as along 
the a-axis. 



compared to the dipole-dipole correlations (|16fl . be- 
cause the neutron scattering cross-section is the product 
of the correlator with the single-dipole contribution, see 
Eq. p7fl in the following Sec. 



8 




-0.05 



0.05 
(q b -Q AF )/(2;i) 



FIG. 6: (Color online). Slices of the normalized correlator Dq 
along the b-direction for three different do-pings. The widths 
of the peaks clearly increases with increasing 



V. INCOMMENSURATE NEUTRON 
SCATTERING PEAKS 

As explained in the introduction, the experimental 
conditions 7 are such as neutrons only interact with the 
projections of electron spins orthogonal to the 6-axis, 
which in terms of the NLSM is just the 7r-field. Using 
Eqs. i|12|) . © and (|16|l . we find the following neutron 
scattering cross section 



oc kt. 



CX 7T 



xfl 



E 



(17) 



where f2 is the area of the sample and Gq is the single- 
dipole contribution given by 



G" = M z 



■. (q ■ e b) 

{q 2 + MD 



1P\ 



(18) 



Experimental results are available for 1%, 1.4% and 1.8% 
doping, shown in Fig. 3 of Ref.lTjand reproduced in Fig.|5J 
Since we assume the presence of a dominant Neel order, 
our results are reliable only at small doping, far from the 
transition to the spin-glass phase. We therefore only con- 
sider x = 0.01 and x = 0.014. It is explained in Ref. Q 
that the small asymmetry experimentally observed for 
1.4% doping [Fig. |Hfb)] arises due to different twinning 
directions (there are four possible equivalent orthorhom- 
bic distortions of the tetragonal lattice) present in the 
crystal. Eq. Ijl7(l clearly shows that the incommensurate 
peak observed in neutron scattering is due to two effects, 
because it is the product of the single-dipole contribution 
Gq with the correlation function Dq. In the limit of non- 
interacting dipoles the correlator is trivial, = 1, and 
hence the maximum of the neutron peak, solely deter- 
mined by the smaller mass, is found at q = M a . However 
in the case of interacting dipoles, there is a second, much 
broader peak in the correlation function which has to 
be taken into account. Fig. [7| clearly shows that neither 
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FIG. 7: (Color online). Normalized neutron scattering cross 
section along the orthorhombic b-axis for x=0.01. According 
to Eq. Gq (a) represents the single dipole contribution 

and Dq (b) the dipole-dipole correlation function. Neither 
of these contributions alone is sufficient to characterize the 
resulting cross-section |Tq| 2 (c) given by their product. For 
comparison with experiments, we convolute |Tq| 2 with a Gaus- 
sian to take into account the finite experimental resolution (d). 
Note that all quantities are normalized to unity. 



of the two contributions alone is sufficient to characterize 
the resulting neutron scattering cross section. Only the 
product l(T7l) gives the right answer. In order to com- 
pare our results with experiments, we take into account 
the finite experimental resolution by a convolution with 
a gaussian of half-width r roso i = 0.005. This convolu- 
tion only leads to a nonzero value at q = but does 
not have any influence on the positions and the widths 
of the peaks, as illustrated in Fig. [7| A comparison with 
experiments, presented in Fig. [SI shows that our simu- 
lations arc in remarkable agreement with experimental 
observations, especially considering the fact that these 
curves do not contain any fitting parameters. Although 
the finite-size effects encountered in our simulations influ- 
ence the correlator Dq (see Sec. lIV|) . they are practically 
absent in the neutron scattering cross section. Simula- 
tions of systems with 140 dipoles (not shown) lead to 
almost identical curves in Fig. |SJ 



VI. DOPING DEPENDENCE OF THE ZERO 
MOMENTUM DM SPIN- WAVE GAP 

Gozar and co-workers 8 have recently measured the DM 
induced spin-wave gap in LSCO at zero momentum using 
Raman scattering. Their result for zero doping A dm = 
D ps 17 ' AcmT 1 « 2.2meV agrees reasonably well with 
the value D ss 2.5meV known from neutron scattering 37 . 
According to Ref. 8, the gap is reduced by 28% to Adm ~ 
12.5cm _1 at 1% doping. As one should expect, the gap 
vanishes at x = 0.02, where the Neel order is completely 
destroyed. Since our approach relies on the existence 
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FIG. 8: Neutron scattering probability for x — 0.01 and 
x — 0.014. The dots correspond to experimental observations 
taken from Fig. 3 in Ref. Q, with normalized intensity. The 
curves represent our simulations convoluted with a gaussian 
to take into account the finite experimental resolution. The 
agreement between simulations (containing no fitting param- 
eters at all) and experiments is remarkable. 



of an ordered state, it does not allow to describe the 
transition to the spin-glass phase at x = 0.02, but for 
x = 0.01 the approximation is well justified. 

Since the anisotropy due to the DM interaction D is 
independent of doping, the relation Adm — D = cM a is 
only valid at zero doping. In order to determine the spin- 
wave gap at finite doping, we consider an excitation with 
zero momentum. In the framework of the NLSM, the 
ground state of the Hamiltonian Q can be represented 
as ft = (n a , n b , n c ) = (sin [7r a (r)] , cos [7r a (r)] , 0), with the 
7?-ficld found in the previous Sees. Since this vector-field 
is small, ir a <C 1, one can expand the trigonometric func- 
tions in powers of 7r a . In order to preserve rotational 
symmetry, we expand only the final answer and use the 
exact expression in our calculations. A zero momentum 
excitation corresponds to a global rotation of the ground 
state n-field by an angle a«l 

n(r) -> n(r, t) = (sin[7r a (r) + a (t)] , cos[7r a (r) + a (t)] , 0) . 



a depends on time, but it is is space-independent, be- 
cause it describes a global rotation. The Lagrangian of 
the system consists of the usual kinetic term, propor- 
tional to (dtn) 2 , and the potential energy E (n) 



L = 



^r^{d t nf-E(n)) 



Jd 2 r{^a 2 -^(M 2 a n 2 a + M 2 n 2 c )} 



(19) 



Without anisotropies, i.e. M a = 0, the system is ro- 
tationally invariant and the potential energy therefore 
independent of the uniform rotation a (Goldstone theo- 
rem). Only the mass terms break the 0(3) symmetry and 
we thus restrict the Lagrangian i|19[l to these terms. The 
anisotropy in c-direction is irrelevant, because n has no 
c-component. Expanding the M a -term in (|19|1 in powers 
of a (omitting zeroth order terms), one finds 



Pi 



Mi 



nzd r 



cos[27r a (r)]d 2 r 



The contribution linear in a disappears, because 
J sm[2ir a (r)]d 2 r = 0. The Eulcr-Lagrange equation of 
motion reads 

n X ±a(t) = -a p s Ml J cos[27r a (r)]d 2 r 

« -a Ps M 2 n(l - 211) , (20) 

where O is the area of the sample, and 

d 2 q ,_ Q ,2 



n 



d z r [n a (r)Y 



(2tt) 2 



d 2 q 



2 g%d^. (2i) 



In order to obtain this last expression, we expanded the 
right hand side of (|20|) in powers of the static field 7r a 
and also used Eq. l(T7|) . The solution of Eq. (|2"U|) gives 
the frequency corresponding to the gap resonance 



*DA1 



c M a (i - n) 



(22) 



In the limit of noninteracting dipoles, = 1, II can be 
calculated analytically. Using Eqs. (fTHfl . and (fT0f> 

one finds II = x^£- (ln(2«/M a ) - 5/4) w 13a;. However, 
the interaction of dipoles is very important. Using the 
correlator DJJ found in Sec. II VI in the numerical integra- 
tion of <J2IJ), wc find that LT ps 25. 5x. At x = 0.01, this 
gives a 26% reduction of the gap, in very good agree- 
ment with the experimental value of 28%. Note that II 
determines the average angle of spin deviation from the 
6-direction, ip 2 = II (see Fig. EJ. Thus at x = 0.01, the 
root mean square value of the angle is ip rms ~ 0.5 rad. 

The reduction of the DM spin-wave gap upon dop- 
ing has recently been considered in Ref. |38l Their ap- 
proach^ is based on the introduction of an effective phe- 
nomenological Lagrangian for the "dipolar field" . Our 
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detailed numerical simulations show that the dynam- with 
ics of the dipoles are mainly diffusive, so it is impossi- 
ble to introduce such an effective Lagrangian. In addi- 
tion, the authors^ suggest a helical spin structure for 
0.02 < x < 0.055. Our results do not support this sce- 
nario. At least in the Neel phase, the spins remain con- 
fined to the afe-plane (apart from a small DM canting) 
and we expect the same for the spin-glass phase. E 



J Bib 



(2) 
Bip 



2V2 



(nx By 



- 2B ■ (D + x ft 



-- I d 2 r (ip f a-riip) (B ■ ft) 



J d 2 r^a{e-p,B -n(S -nj}ip , (24) 



VII. UNIFORM MAGNETIC 
SUSCEPTIBILITIES 



The magnetic susceptibility of undoped LCO, includ- 
ing its temperature dependence is well understood 3 - 
(see also a recent work by Silva Neto et alS^). There 
are four mechanisms that contribute to the magnetic 
susceptibility of LCO: 1) isotropic atomic core dia- 
magnetism, 2) anisotropic van Vleck paramagnetism, 
3) anisotropic quantum Hciscnbcrg model paramag- 
netism, and 4) anisotropic paramagnetism related to the 
relativistic DM interaction. In the present work, we con- 
sider the doping dependence of the uniform susceptibili- 
ties at zero temperature. Due to the anisotropies, there 
are three different susceptibilities, Xa, Xb, and Xc, corre- 
sponding to directions of the magnetic field along the a-, 
b-, and c-axis, respectively. 

The magnetic susceptibilities of lightly doped LSCO 
have been measured by Lavrov et al~. Fig. 1 of Ref. @ 
clearly shows that their doping dependence is strongly 
anisotropic: Xa and Xc change only slightly with dop- 
ingifi, whereas Xb varies from about 1.7 • 10~ 7 emu/g at 
x = to 4 • 10 -7 emu/g at x = 0.01. In comparison, the 
perpendicular susceptibility of the Heisenberg model is 
equal to 



X± « 0.5/(87) 



16J 



1.6 • 10 emu/g 



Here we have restored the gyromagnetic ratio, g s w 2, 
and the Bohr magneton [Ib (throughout this paper, we 
use units where g s fJ-B = 1) and substituted the real den- 
sity of the compound. The variation of the longitudinal 
susceptibility Xb (it is longitudinal because the field is 
directed along the Neel direction) at a tiny 1% doping is 
comparable or even slightly larger than x±- Similar to 
the neutron scattering considered in the previous Sees., 
this is an indication of a noncollinear spin structure. 

In the absence of a magnetic field, the Hamiltonian of 
the system is given by Eqs. ©, (0J, and JSJ. The addi- 
tional terms describing the interaction with the magnetic 
field read 



where {, } stands for the anticommutator, p = — iV and 
we set g s fJ.B = 1- Note that the field n and the spinor 
ip are functions of r. The first term in Esn is the usual 
magnetic interaction in the NLSM^S^i and the second 
term is due to the DM induced weak ferromagnetisrn^i* 3 ^. 

Eg^ describes the interaction between the hole and 
the component of the magnetic field parallel to the lo- 
cal direction of the n-field. The physical origin of this 
term is very simple. In the Heisenberg model (undoped 
case), the projection of the total spin along ft is zero, be- 
cause electrons with spin "up" are compensated by elec- 
trons with spin "down" . For nonzero doping, "up" and 
"down" -spins no longer compensate and therefore lead 
to a nonzero spin-projection along ft. This spin interacts 
in the usual (Zeeman) way with the magnetic field. For 
example, ■ ri^ = +1 implies that the hole is created 

on the f-sublattice and hence AS n = +1. 

(2) 

The origin of E B ^ is more delicate. This term describes 
the interaction between the hole and the component of 
the magnetic field orthogonal to the local direction of 
ft. The physical origin of this interaction is the follow- 
ing: the locally transverse magnetic field tilts the spins 
in the antiferromagnetic background, which allows the 

holes to hop more easily and therefore leads to a gain in 

(2) 

kinetic energy. This is why the -E^-term contains the 
momentum, p = — zV, and is proportional to the hopping 

(2) 

matrix element g — Zt. To derive the coefficient in E B ^ 7 
one has to calculate the gain in the kinetic energy and 
then expand it for a hole localized in a particular pocket, 
q = (±7r/2, ±7r/2) + p. This why the unit vector e, or- 
thogonal to the corresponding face of the MBZ appears 
in this expression. However, because in the problems we 

consider, momenta are always small, p< 1, we can safely 
(2) 

neglect the I^A-term. We nevertheless present the term 
in Eqs. lll'-il) and (|24fl for the sake of completeness and 
for future studies. 

To the best of our knowledge, the interaction between 
the hole and the magnetic field has never been written 
in the form l|24() before. It is therefore useful to convince 
oneself of its validity by recalculating a well known re- 
sult. A nice example is the calculation of the perpendic- 
ular susceptibility of a single immobile hole in a quantum 
antifcrromagne^Si^. Since the hole is located on a par- 
ticular site of the lattice, ip^a ■ ni/) — ±5 (r), and because 
9 = (no hopping), the total energy reads 



E f 



E 



(i) 

Bib 



E 



(2) 

Bib 



(23) 
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Here the terms in square brackets are due to Eq. J7J, 
and the last term is due to Ehl in l|24|) . Variation of the 



energy with respect to yields 



B 



(-A + M 2 ) 

which in momentum representation has the solution 



'(r) = ±-^(r) 
2p s 



B 



1 



M 2 



Substituting this solution in the total energy and using 
the ultraviolet regularization r > 1 lattice spacing, i.e. 
q < 1, then immediately yields the "impurity" suscepti- 
bility 



Xn 



87rp s \M 



which agrees with previous results obtained by different 
methods^. 

Let us now turn to the realistic situation of LSCO, 
where mobile holes are trapped by Sr ions ( "impurities" ) . 
To calculate the susceptibility, it is sufficient to consider 
a single "impurity" and then multiply the result by the 
concentration x, because the interaction of "impurities" 
considered in Sees. IIIII and IIVI gives corrections of second 
and higher order in doping, which we cannot calculate 
without uncontrolled approximations, see also comment 
below Eq. (|28[1 . In what follows, we set M a = every- 
where except in the logarithmically divergent integrals, 
where we use it as an infrared cutoff. Accordingly, the 
term involving D + in Eq. I|24|) is also ignored. 



A. Longitudinal magnetic field, Xb 

In the case of a longitudinal magnetic field B — Be b , 
the "impurity" energy is given by Eqs. (7J|, l|23|) . and i|24|) 



\ I d 2 r{ ,,, 



X±B 2 



+ 2V2gp(r) 

-I 2 



V6»(r 



(e 6 • V) 9 (r) 



■0) 



— B 



dH+d x -e(r)]p(r) , (25) 



where we represent the vector d — ^<jty as d — d" e b +d A 
Variation of the energy with respect to 9 yields 



-A0(r) 



M ( e b x d x ) (e b ■ V) + 



Bd A 
~2p~s 



p(r) = . 



In momentum representation, the solution of this equa- 
tion reads 



iM (e b x d^ e b ■ q + 



Bd A 
~2p~ s 



Pq_ 

«2 



(26) 



where the first term, which is independent of the mag- 
netic field B, coincides with 0. Note that analogous 
to neglecting the mass terms, we also omit the \^B 2 9 2 
contribution in the solution l|26|l . because it only leads to 
a small modification of the infrared cut-off. Substitution 
of this solution in the total energy (|25|l yields 



E = -{d A 



B 2 
4irp s 



XJ-9 



In 



4irp s 



0.6k 



(27) 



with the kinematical constraint (d") 2 + (d ± ) 2 = 1. The 
energy is minimized for a non-zero d" , i.e., a longitudinal 
magnetic field leads to a squeezing of fa = e b x d. The 
minimum is obtained for 



Kp s 

i i 



B 



Substitution of this expression in l|27l) gives the "impu- 
rity" energy and allows to find the "impurity" suscepti- 
bility. After multiplication by the hole concentration x, 
we find the following variation of the bulk susceptibility 



( c 2 



9 



\8irp 2 
100x X ± 



2TT P i 



In 



0.6k 



7TC 



2g 2 K 2 



XXA- 



(28) 



The last term in the above expression is the most impor- 
tant one. Depending on the value of k, it gives 80% - 
90% of the total result. From data on hopping conduc- 



tivity 



.32.33 



we know that 0.2 < k < 0.3. Evaluation of 



Eq. l(2^|) for k = 0.3 gives 5\b ~ 75xx_i_ and with k = 0.2, 
we find Sxb ~ 130.txx. Similar to the uncertainty in k, 
there is also some variability in the coupling constant g. 
According to Refs. IlllllSl we take g ss 1, but we believe 
that g w 1±0.2 is quite possible. Despite these uncertain- 
ties, our result unambiguously shows that there is a huge 
doping dependent variation of the longitudinal suscepti- 
bility due to the noncollinear spin structure of the spi- 
ral. Our estimation of 8\b is in excellent agreement with 
experimental datai We would like to emphasize that 
Eq. (|28|l has been derived in the fully controlled linear 
in x approximation. Unfortunately, higher order doping 
terms cannot be calculated without uncontrolled approx- 
imations. Eq. (|28|l is therefore justified when x is well 
below the transition to the spin-glass phase at x = 0.02. 
Since our results on neutron and Raman scattering are 
also derived in the linear in x approximation, and we 
know from comparison with experiments that they are 
valid up to x ~ 0.015, we expect Eqs. (P%)l and (|3*T1) to be 
valid in this region as well. In contrast to this huge mod- 
ification of the longitudinal susceptibility upon doping, 
practically no variation of the transverse susceptibilities 
has been observed experimentally 9 . 



12 



B. Transverse magnetic field, \a and \ c 

For a transverse magnetic field B _L <={,, the "impurity" 
energy given by Eqs. J7J), 1(25^. and (J23J reads 



^ = 5/ [ V ^ (r) ] 2 

+ 2V2.9 (e 6 x <f) p (r) (e 6 • V) (r) 



B-9{r) - B ■ 9 (r) d x • 6 (r) p (r) I . (29) 



In this case, c?" = and d = d^, which maximizes the 
"impurity" dipole moment and hence minimizes the en- 
ergy. Let us write the 6*-field as 



6 (r) = 6f(°) (r)+0"W ( r ) , 



(i) 



(30) 



where 9 ^ is the solution for B = given by Eq. © and 
9 W is the perturbation induced by the magnetic field. 
Performing the variation of the energy (|29() with respect 
to 6 W yields the equation 



A(9 (1) (r 
2p s I 



^ ■ 6» (r) 



(r)]}p(r) 



. 



Using the explicit form of 9 ^ given by Eq. © (we also 
set M a = 0) we then find the magnetic field induced part 
of the spiral 



9(1) 



n 2 M 

■1 

2TTp s 

2 k 



Q 

2nq 2 



4k 



y/4:K 2 + q 2 (16k 2 + 9 2 ) 3/2 y/WK 2 + q 2 



where we have used the fact that 9 cx e a . The term 

(^B ■ 9^j in (|29|l is thus only present for a magnetic field 

directed along the a-axis. Such a small variation of the 
susceptibilities upon doping is quite consistent with ex- 
periments^. 

VIII. CONCLUSION 



In the present work based on the spiral solution of the 
extended t-J model, we explained the following proper- 
ties of underdoped La2-a:Sr K Cu04. 

1) The pinning of the incommensurate magnetic struc- 
ture to the orthorhombic 6-direction observed in neutron 
scattering in the insulating phase, x < 0.055. The pin- 
ning is due to the anisotropy of the diagonal hopping 
matrix element t' , see Sec. ITT1 

2) The positions and shapes of the incommensurate elas- 
tic neutron scattering peaks in the Neel phase, x < 0.02. 
Experimental data are presented in Fig. [SJ together with 
our theoretical curves, containing no fitting parameters 
at all. The agreement between theory and experiments 
is quite remarkable. 

3) The doping dependence of the Dzyaloshinski-Moriya 
induced spin- wave gap in the Neel phase, see Eq. J22J). 
According to our calculation, at 1% doping, the gap is 
reduced by 26% compared to its value in the undoped 
compound. This is in very good agreement with the ex- 
perimentally observed reduction of 28%. 

4) The doping dependence of the uniform magnetic sus- 
ceptibilities at zero temperature, see Eqs. (|28f) and 13111 . 
This explains the tremendous variation of the longitudi- 
nal magnetic susceptibility \b and the very weak change 
in the transverse susceptibilities Xa and Xc- 



Substitution of this solution together with @ in Eq. I|30l) 
and then in Eq. (|29|l gives the total energy and hence 
the "impurity" susceptibility. After multiplication by the 
hole concentration x, we find the variation of the bulk 
susceptibilities 



5Xa = 



^ 2 lnf^ + 0.0043.^ 2 



An 
Wxx± , 



M a 



4:Tr 3 p s 



k 2 M 2 

dxc = 0.0043 — ^— x ~ 0.02x X ± 



(31) 
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